METHOD FOR MEASURING BONE MINERAL DENSITY BY USING X-RAY IMAGE 

Field of the Invention 

5 The present invention relates to a method for measuring a bone 

mineral density using an X-ray image. More particularly, the present 
invention relates to a method for measuring a bone mineral density 
by removing an X-ray absorption effect due to soft tissues contained 
in a bone portion of an X-ray image, when the bone mineral density 
Z 10 is measured using the X-ray image, and a storage medium being readable 
jTj! on a computer storing a program to implement the method. 

i 

" Description of the Prior Art 

J=J 15 Osteoporosis is a wide spread medical condition that affects 

f ; the middle-aged and older populations. Especially, the condition 
A is prevalent in postmenopausal women. Osteoporosis is characterized 

by an abnormal loss in bone mineral content, which leads to a tendency 
toward non-traumatic bone fractures and to structural deformations 
20 of bones. However, effective therapy for osteoporosis has not been 
developed yet. Accordingly, it is important that a method for easily 
and inexpensively diagnosing the osteoporosis should be developed 
for the prevention of deterioration of osteoporosis and early stage 
treatment of osteoporosis. 
25 Bone mineral density is one of important factors for diagnosing 

osteoporosis. Various bone mineral density measurement methods have 
been developed. 
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Quantitative computed tomography (QCT) provides a 
three-dimensional bone density image and thus provides separate 
estimations of cortical and trabecular bone densities. Based on 
the three-dimensional bone density image, QCT method can provides 
5 a structural stability of a bone to some extent. However, there 
are some limitations to use the QCT as a routine screening tool 
for osteoporosis because the price of the QCT equipment is very 
high and the radiation dose of a QCT scan is generally several hundred 
times larger than that of a plain x-ray imaging. 
|=J 10 The most widely used method for measuring bone mineral density 

f and for follow-up study of osteporosis patients is a dual energy 
;~; x-ray absorptiometry (DEXA). The precision error of the DEXA in 

determining bone mineral density is reported to be about few percents . 
Furthermore, the radiation dose of a DEXA scan is very small compared 
^ 15 with a QCT scan. 

;j; Broadband ultrasound absorption ( BUA) is also used for measuring 

: * bone mineral density. However, BUA is not so accurate in determining 

bone mineral density compared with other equipments. Nevertheless 
the validity of BUA for osteoporosis study does not diminish because 
20 recent studies have reported that BUA is somewhat relevant to the 
mechanical strength of bone. 

In spite of the various methods prescribed above, such as QCT, 
DEXA, and BUA, methods for measuring bone mineral density and 
diagnosing osteoporosis by using an x-ray image are developed steadily 
25 in the practical point of view. The reason is that most hospitals 
are equipped with an x-ray imaging system basically, so there is 
no need for an extra cost to purchase a new bone mineral density 
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measuring system. Moreover, the high quality of an x-ray image enables 
trabecular pattern analysis. Particularly, it is believed that the 
trabecular pattern contains useful information about the fracture 
risk due to osteoporosis. In this sense, many researchers study 
trabecular patterns to extract useful information related to the 
fracture risk. 

Up to now, various methods for measuring bone mineral density 
by using an x-ray image have been developed. However, bone image 
in the X-ray image contain significant amount of x-ray absorption 
effects due to overlapping soft tissues. Actually, in the 
anterior-posterior x-ray image of a femoral neck, the amount of 
x-ray absorption due to overlapping soft tissues is comparable to 
the amount of x-ray absorption due to the femoral neck itself. By 
the way, x-ray absorption due to overlapping soft tissues can not 
be easily removed in the simple x-ray imaging scheme, which leads 
to a relatively low accuracy in the bone mineral density measured 
by using an x-ray image. 

Summary of the Invention 

It is, therefore, an object of the present invention to provide 
a method for measuring a bone mineral density by removing an X-ray 
absorption effect due to soft tissues contained in a bone portion 
of an x-ray image, when the bone mineral density is measured using 
the X-ray image, and a storage medium being readable on a computer 
storing a program to implement the method. 

In accordance with an aspect of the present invention, there 



is provided a method for measuring a bone mineral density, by use 
of an x-ray image, in a bone mineral density measuring system, 
comprising the steps of: (a) obtaining an X-ray image of a bone; 
(b) setting a region of interest on the obtained X-ray image of 
5 bone; (c) calculating a background trend due to soft tissues, at 
a bone portion within the set region of interest; and (d) calculating 
an index of the bone mineral density by removing the background 
trend due to the soft tissues, at the bone portion within the set 
region of interest. 

10 In accordance with another aspect of the present invention, 

there is provided a storage medium being readable on a computer 
storing a program to implement a method for measuring a bone mineral 
density, by use of an x-ray image, in a bone mineral density measuring 
system, comprising the steps of: (a) obtaining an X-ray image of 

15 a bone; (b) setting a region of interest on the obtained X-ray image 
of bone; (c) calculating a background trend due to soft tissues, 
at a bone portion within the set region of interest ; and ( d ) calculating 
an index of the bone mineral density by removing the background 
trend due to the soft tissues, at the bone portion within the set 

20 region of interest. 

Brief Description of the Drawings 

The above and other objects and features of the present invention 
25 will become apparent from the following description of preferred 
embodiments taken in conjunction with the accompanying drawings, 
in which: 
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Fig. 1 shows a bone mineral density measuring system to which 
the present invention is applied; 

Fig. 2A is a flowchart illustrating a method for measuring 
a bone mineral density using an X-ray image in accordance with an 
5 embodiment of the present invention; 

Fig. 2B is a detailed exemplary flowchart describing a step 
for obtaining the X-ray image shown in Fig. 2A; 

Fig. 2C is a detailed exemplary flowchart describing a step 
for calculating a background trend due to soft tissues within a 
10 radius portion in a region of interest shown in Fig. 2A; 

Fig. 3 is an exemplary view showing an aluminum step-wedge 
used to calibrate the X-ray image in accordance with the present 
invention; 

Fig. 4 is an exemplary view of the X-ray image obtained after 
15 imaging a step-wedge and a wrist together in accordance with the 

present invention; 

Fig. 5 illustrates an exemplary relationship between the 

thickness of the step-wedge and the average gray-level of the 

step-wedge shown in Fig. 4; 
20 Fig. 6A illustrates the X-ray image calibrated by the step-wedge 

and a rectangular region of interest selected to measure the bone 

mineral density of a radius bone in accordance with the present 

invention; 

Fig. 6B illustrates a wrist cross-sectional view taken along 
25 a horizontal line 1 of the region of interest; 

Fig. 7A illustrates a gray-level profile taken along the 
horizontal line 1 shown in Fig. 6A; and 
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Fig. 7B illustrates a gray-level profile along the horizontal 
line 1 after the background trend due to the overlapping soft tissues 
in the radius portion is removed. 

5 Detailed Description of the Invention 

Referring to Fig. 1, there is shown an exemplary diagram 
illustrating a bone mineral density measuring system for measuring 
a bone mineral density in accordance with the present invention. 

10 As shown , the system for measuring a bone mineral density includes 

an input /output device 11 for inputting or outputting data required 
for measuring the bone mineral density by a user, main/auxiliary 
storage devices 12 and 13 for storing various data necessary for 
a process of measuring the bone mineral density by using an x-ray 

15 image, and a microprocessor 14 for controlling the input/output 
device 11 and the main/auxiliary storage devices 12 and 13 and 
controlling full operations for measuring the bone mineral density 
by using an x-ray image. 

The input /output device 11 includes a monitor, a printer, and 

20 an x-ray film scanner for digitizing an x-ray film. If a digital 
image sensor (CCD or CMOS sensor) is used instead of an x-ray film, 
the digital image sensor is included in the input/output device 
11 . 

Through the above system, the method for measuring the bone 
25 mineral density by using the x-ray image is carried out in accordance 
with the following procedures. At first, the x-ray image is inputted 
into the input/output device 11 with a program including a process 
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of the following Figs. 2A to 2C incorporated in the microprocessor 
14. When the program is executed, the bone mineral density is measured 
by using the x-ray image according to the program. 

The process is described in greater detail hereinafter with 
5 reference to Figs. 2A to 7. 

Fig. 2A is a flowchart illustrating a method for measuring 
bone mineral density by using an x-ray image in accordance with 
the present invention. Fig. 2B is a detailed exemplary flowchart 
showing a process for obtaining the x-ray image according to the 
; ; f 10 present invention. Fig. 2C is a detailed exemplary flowchart showing 
a process for calculating a background trend due to soft tissues 
;^ within a radius portion in a region of interest. 

H Fig. 3 is an exemplary view illustrating an aluminum step-wedge, 

which is utilized to calibrate the x-ray image, in accordance with 
[II 15 the present invention. The aluminum step-wedge is utilized to adjust 
! s * the variation of the image characteristics dependent on the x-ray 
!'=& generation unit, the types of a screen or a film, the development 
condition of the film, and the characteristics of an x-ray film 
digitizer (the characteristics of the CCD or CMOS sensor, in case 
20 of a digital image sensor) . 

Fig. 4 is an exemplary view of the X-ray image obtained after 
imaging a step-wedge and a wrist together in accordance with the 
present invention. Fig. 5 illustrates an exemplary relationship 
between the thickness of the step-wedge and the average gray-level 
25 of the step-wedge shown in Fig. 4. Fig. 6A illustrates the X-ray 
image calibrated by the step-wedge and a rectangular region of interest 
selected to measure the bone mineral density of a radius bone in 
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accordance with the present invention. Fig. 6B illustrates a wrist 
cross-sectional view taken along a horizontal line 1 of the region 
of interest. Fig. 7A illustrates a gray-level profile taken along 
the horizontal line 1 shown in Fig. 6A. Fig. 7B illustrates a 
5 gray-level profile along the horizontal line 1 after the background 
trend due to the soft tissues is removed. 

At first, as shown in Fig. 2A, the x-ray image of a wrist is 
obtained to measure the bone mineral density in accordance with 
the present invention (step 21). When the x-ray image is obtained, 

10 the conditions for obtaining the x-ray image shall be kept constant 
as possible. In particular, a tube voltage (kVp) and a tube current 
(mA) of an x-ray generation unit should be fixed to specific values. 
For example, when the x-ray image of the wrist described above 
is obtained, the tube voltage and the tube current of the x-ray 

15 generation unit shall be 50kVp and 50mA, respectively. The x-ray 
image is obtained by digitalizing an x-ray film, which was obtained 
by a simple x-ray imaging process, using an x-ray film scanner. 
When a digital image sensor is used, the x-ray image is obtained 
directly without a step of using the film scanner. When the x-ray 

20 image is obtained, its spatial resolution shall be 200PPI (pixels 
per inch) and each pixel shall be of an 8-bit gray level. 

Fig. 4 is an exemplary view of an x-ray image of a wrist, which 
is obtained by digitalizing an x-ray film. Each pixel of the x-ray 
image consists of a gray-level associated with the amount of x-ray 

25 absorption by the body. The amount of x-ray absorption is determined 
by a density or thickness of bone and soft tissues. Generally, the 
x-ray image of bone also includes x-ray absorption effects due to 



overlapping soft tissues. Accordingly, to accurately measure bone 
mineral density using the x-ray image, it is necessary to remove 
x-ray absorption effects due to the overlapping soft tissues from 
the x-ray image. 

5 In the meantime, the characteristics of the x-ray image, such 

as brightness, contrast, etc. are sensitive to the change of x-ray 
imaging conditions . In other words , the brightness , i.e. gray-level, 
of the x-ray image at each pixel cannot be an absolute measure for 
the amount of x-ray absorption. Therefore, when the bone mineral 
i; 10 density is measured by using the x-ray image, it is necessary to 

3 calibrate a variation of image characteristics due to the change 

y 

^ of the x-ray imaging conditions and to quantify the amount of x-ray 
*jj absorption at each pixel. 

In order to calibrate the variation of characteristics of the 
y 15 x-ray image due to x-ray imaging conditions, and quantify the amount 
a of x-ray absorption at each pixel, a metallic step-wedge composed 
=& of several steps is used as shown in Fig. 3. Aluminum or copper 
is suitable for a material of the step-wedge. Generally, the material 
and thickness of the step-wedge are properly selected according 
20 to a body portion where the bone mineral density is measured. 
Conventionally, the material and thickness are selected in such 
a way that the amount of x-ray absorption by the highest step of 
the step-wedge is slightly larger than the maximum amount of x-ray 
absorption by the body portion where bone mineral density is to 
25 be measured. When the bone mineral density is measured at the wrist, 
the appropriate maximum thickness of the aluminum step-wedge is 
12mm. For a step-wedge containing heavy metals such as copper, etc . , 
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the maximum thickness shall be reduced in consideration that the 
x-ray absorption becomes larger. The step-wedge shown in Fig. 3 
is made of an aluminum plate of which a bottom surface is a square 
with a dimension of 40mm x 40mm and thickness is 12mm, processed 

5 into radial steps having eight steps. In the aluminum step-wedge 
described above, the heights of steps are 1.5, 3.0, 4.5, 6.0, 7.5, 
9.0, 10.5 and 12.0mm, in an increasing order. 

As shown in Fig. 2B, in order to calibrate a characteristic 
variation of the x-ray image due to the x-ray imaging conditions, 

10 and to quantify the amount of x-ray absorption at each pixel, an 

i-S x-ray image of the radial step-wedge and a wrist is obtained (step 

i.r= 

211). Fig. 4 is an exemplary view of the x-ray image where the 

□ 

step-wedge and the wrist are imaged together. In Fig. 5, 
t two-dimensional data between the thickness of the step-wedge and 

? 15 the average gray-level of the step-wedge shown in Fig. 4 are denoted 

by solid squares . In Fig . 5 , the average gray-level , when the thickness 

Q 

■ A of the step-wedge is zero, corresponds to the average of background 
gray-levels around the step-wedge. 

Then, as shown in Fig. 2B, the full x-ray image is calibrated 
20 by using above two-dimensional data between the thickness of the 
step-wedge and the average gray-level of the step-wedge (step 212). 

In order to calibrate the full x-ray image by using the two-dimensional 
data, average gray-level at any thickness of the step-wedge should 
be calculated. In order to calculate the average gray-level at any 
25 thickness of the step-wedge, the two-dimensional data shown in Fig. 
5 are fitted to an appropriate function. In Fig. 5, the gradient 
of the average gray-level increases gradually at first and then 
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decreases gradually as the thickness of the step-wedge increases. 
One of the functions with the characteristics described above is 
a tangent hyperbolic function. Thus a fitting function of "f(t) 
= a+b*tanh (c*t+d)" is used, where "f (t)" represents the gray-level, 
5 "t" represents the thickness (mm) of the step-wedge, and "a", "b", "c" 
and "d" represent the fitting parameters, respectively. Here, the 
function "f(t)" is symmetric with respect to a point of "t = -d/c". 
However, as shown in Fig. 5, the gray level profile of the step-wedge 
is not completely symmetric. Thus a fitting of data in Fig. 5 to 

j.j, 

10 a single fitting function "f(t)" may create considerable fitting 

■■■-i 

errors. Accordingly, in the present invention, the data are divided 
'|f into two parts and then data of each part are fitted to the function 
'f } "f(t)". The fitting function "f(t)" has four fitting parameters, 
so there needs to have more than four data in each fitting part. 
s ff 15 The first part consists of six data with heights of 0mm to 7.5mm. 

The second part consists of six data with heights of 4 . 5mm to 12 . 0mm. 
! * Here, the overlap of the two parts is to softly couple fitting 

results. The fitting is carried out in each part. The 
Levenberg-Marquardt fitting method is used as a fitting method in 
20 the present invention. 

Here, "fi(t>" denotes the fitting result in the first part, 
and "f 2 (t)" denotes the fitting result in the second part, and the 
results are shown in Fig. 5 as a solid line and a dotted line, 
respectively. Then, the two fitting functions are coupled together 
25 to provide a final fitting function "F(t)" as follows. At first, 
"F(t) = fi(t)" for "t <= 4.5", "F(t) =f 2 (t)" for "t >= 7.5", and 
"F(t) = x*fi(t) + (l-x)f 2 (t)" for "4.5 < t <7.5". Here x = (7.5-t)/3. 
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The function "F(t)" is a simply increasing function, so the inverse 
function F _1 (g) is well defined. Here "g" denotes the gray-level. 
Now, the gray-level "g" of each pixel of the x-ray image is calibrated 
as following procedures by using the final fitting function "F(t)". 
In case of "g >= F(12)", the calibrated gray-level is "255". In 
case of "g <= F(0)", the calibrated gray-level is 0. In other case, 
the calibrated gray-level is an integer part of "256*F -1 (g) /12" . 

Fig. 6A shows a calibrated image of Fig. 4, which was obtained 
via the calibration method described above. 

Then, as shown in Fig. 6A, in order to remove the X-ray absorption 
effect due to soft tissues contained in the radius portion of the 
X-ray image, a rectangular region of interest is set on the radius 
portion ( step 22 ) . A size of the region of interest is 250x150 pixels, 
and the region of interest should contain soft tissue regions to 
the left and right sides of the radius. 

Then, the background trend due to soft tissues is calculated 
within the radius portion in the region of interest (step 23). 

In order to help the understanding of the method for calculating 
the background trend due to soft tissues, a wrist cross-section 
"A" along a crossing line 1 of Fig. 6A is shown schematically in 
Fig. 6B. An inner dark portion in the cross-section "A" is a radius 
portion and an outer bright portion is a soft tissue portion . Generally , 
the cross-section of a long bone such as a radius bone is approximately 
a disc shape. Accordingly, the cross-section of the radius bone 
is simplified as a disc shape hereinafter. As shown in Fig. 6B, 
the cross-section "A" is decomposed as follows: (cross-section "A") 
= (cross-section "B" ) + (cross-section "C") - (cross-section "D"), 
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where the cross-section "B" is equal to the cross-section "A" except 
that the radius bone is replaced with soft tissues, the cross-section 
"C" is equal to the bone portion in the cross-section "A", and the 
cross-section "D" is equal to the cross-section "C" except that 

5 the radius bone is replaced with soft tissues. On the other hand, 
the gray-level profile along the line 1 in Fig. 6A is drawn by a 
solid line in Fig. 7A. In Fig. 7A, the horizontal axis denotes a 
coordinate of pixel and the vertical axis denotes a gray-level. 
In Fig. 7A, intervals "ai-bi", "bi-Ci" and "ci-di" are a soft tissue 

10 region, a radius bone region, and another soft tissue region, 
respectively. 

In the interval "bi-Ci" of Fig. 7A, it is impossible to exactly 
calculate the background gray-level trend due to soft tissues from 
the simple x-ray imaging scheme. Accordingly, an approximate method 

15 is introduced to calculate the background trend. In the present 
invention, we use an interpolation method where gray-level profiles 
of the two soft tissue regions "ai-bi" and "ci -dj." are interpolated 
into the radius region "bi-Ci" to set it to a background trend. 
For the sake of interpolation, a dif f erentiable fitting function 

20 is selected (step 231) . Generally, a polynomial function is suitable 
for the fit. However, in the present invention, only a 4 th order 
polynomial is used for the fit. The interpolation is accomplished 
by fitting data of the soft tissue regions "ai-bi" and "Ci-di" into 
the 4 th order fitting function through a Levenberg-Marquardt fitting 

25 method. The interpolation result in the radius bone region is set 
to the background trend due to soft tissues. That is, the gray-level 
profiles of the soft tissue regions adjacent to the radius region 
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are interpolated into the radius region using the chosen fitting 
function to set the background trend (step 232) . Seeing this method 
in detail, the method includes a step of obtaining the gray-level 
profile across the radius bone in the region of interest, a step 

5 of dividing the gray-level profile into the bone portion and the 
soft tissue portions, and a step of interpolating the gray-level 
profiles of the soft tissue portions into the bone portion with 
the fitting function to set the fitting result to the background 
trend due to the soft tissues in the bone region. "LI" of Fig. 7A 

10 is the background trend interpolated via the above method. This 
background trend corresponds to the gray-level profile due to the 
cross-section "B" of the Fig. 6B. Now, the background trend is 
calculated with the cross-sectional line 1 moved to all the rows. 
If the background trend has been set, the background trend 

15 due to the soft tissues is removed from gray levels of the radius 
region to calculate an index for the bone mineral density (step 
24). The step 24 includes a step of removing the background trend 
in the gray level of each pixel of the bone region, a step of calculating 
an average gray level (<G>) in the bone region after removing the 

20 background trend, a step of calculating a weighted average (P) of 
radius bone widths within the region of interest, and a step of 
establishing an index for the bone mineral density as <G> + c 0 P. 
Here, c 0 is an unknown constant that should be determined. 

Fig. 7B shows the gray-level profile after removing the 

25 background trend in the radius region "bi-Ci" . The gray-level profile 
of Fig. 7B may be approximately the same as the gray-level profile 
due to the cross-section "C" in Fig. 6B after removing the gray-level 
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profile due to the cross-section "D". The bone mineral content of 
the radius bone is determined by a sum of the gray-level profiles 
due to the cross-section "C" in Fig. 6B. Thus, in the radius portion 
"bi-ci", the sum (G;l {C) ) of the gray-level profile due to the 
5 cross-section "C" is obtained by adding a sum of the gray-level 
profile due to the cross-section "D" to a sum of the gray-level 
profile of Fig. 7B. 



In Eq. 1, "G in " denotes the gray-level profile in Fig. 7B, "G in (D) " 
denotes the gray-level profile due to the cross-section "D", and 



^ gray-level profile due to the cross-section "D" should be known 

\* 15 and in order to know this gray-level profile the shape of the 
!'* cross-section "D" should be known. However, it is impossible to 

know an exact shape of the cross-section "D" from the X-ray image. 

Accordingly, the cross-section "D" of the radius bone is approximated 

to a disc based on the fact that it is anatomically close to the 
20 disc and a diameter of the disc is defined as pi = |ci - bi| . Then, 

"Gi (C) " is approximately equal to Eq. 2. 




Eq. 1 



;;f! io 
■I 



'n" denotes an index for pixel. In order to calculate "Gm 1 



the 




Eq. 2 
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In Eq. 2, "c 0 " is a constant. 
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Now, after calculating Eq. 2 with respect to all the 
cross-sectional lines within the region of interest, the sura of 
all the Gi (C) is given by the following Eq. 3. 



2 



Eq. 3 



The left-hand side of Eq. 3 is an index for total radius bone 
mineral content within the region of interest. 

Now, the left and right-hand sides of Eq. 3 are divided by 
the area of the radius bone A = ^p, within the region of interest, 

and the obtained value is determined as the index for the radius 
bone mineral density (G (C) ). 



In order to apply Eq. 4 to the measurement of the radius bone 
mineral density, the unknown constant "c 0 " shall be determined at 
first. In order to determine the constant "c 0 ", a large number of 
wrist x-ray images are obtained from different female subjects. 
After setting a region of interest for each wrist x-ray image, the 

first term <G>= — YYG ln and the second term P = — Vpf in the 



right-hand side of Eq. 4 are calculated. The second term "P" can 
be replaced with a formof Yp, •fo/Ypk). Then "P" represents a weighted 




Eq. 4 
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average of the radius bone widths within the region of interest 
described previously. Generally, the weighted average is not so 
different from a simple average, so the "P" can be chosen as a simple 
average of "pi" within the region of interest, i.e. the average width 

5 of the radius bone within the region of interest. 

Each female subject has undergone the bone mineral density 
(BMD) measurement within the same radius portion by using a DEXA 
equipment. Then, the constant c 0 is set to a value that minimizes 
the least-squares fit error between BMD and "<G> + c 0 P" . As for 200 

10 female subjects, the constant c 0 set by the least-squares fit is 
0.10. 

Now, as the necessary constant c 0 is determined, the index for 
the radius bone mineral density is given by <G> + c 0 P. Generally, 
the value <G> is several times larger than the value c 0 P, and the 

15 average subject-to-subject variation of the value P is about 4%. 
On the contrary, the average subject-to-subject variation of the 
value <G> is about 40%. Accordingly, the index for the bone mineral 
density can be approximately set by only <G>. 

As described in the concrete embodiment, the numbers or images 

20 used in the embodiment can be replaced with other ones to improve 
the performance of the method by the present invention. The key 
point of the present invention is to measure the bone mineral density 
by removing the x-ray absorption effect due to overlapping soft 
tissues in the bone image, on measuring the bone mineral density 

25 by using the x-ray images. The prior art to measure bone mineral 
density by using x-ray images didn't remove the x-ray absorption 
effect due to the overlapping soft tissues on the bone image, so 
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there remained many possibilities of resulting in a lot of errors. 
However, in the present invention, a gray-level profile of the 
soft tissues surrounding the bone portion is interpolated into the 
bone portion. Then, in the bone region, the interpolated gray-level 
5 profile is removed from the gray-level profile of the bone portion. 
Furthermore, in the present invention, the method of modeling the 
cross-section of the bone to a disc shape is used to compensate 
for the error due to the size or shape variation of the bone. 

The method of the present invention can be implemented by a 
10 program and stored in a computer-readable medium, e.g., a CD-ROM, 
a RAM, a floppy disk, a hard disk, an optical magnetic disk, etc. 

Although the preferred embodiments of the invention have been 
disclosed for illustrative purposes, those skilled in the art will 
appreciate that various modifications, additions and substitutions 
15 are possible, without departing from the scope and spirit of the 
invention as disclosed in the accompanying claims. 
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